Stochastic Gravitational Wave Production After Inflation 



O 
O 



(N 
(N 

(N 
> 

i> : 

\o • 

o : 

O ■ 

^ ' 

Or 

6 ■ 
a: 



13 



Richard Easther and Eugene A. Lim 
Department of Physics, Yale University, New Haven CT 06520, USA 
Email: richard.easther@yale.edu eugene.lim@yale.edu 

In many models of inflation, the period of accelerated expansion ends with preheating, a highly 
non-thermal phase of evolution during which the inflaton pumps energy into a specific set of mo- 
mentum modes of field(s) to which it is coupled. This necessarily induces large, transient density 
inhomogeneities which can source a significant spectrum of gravitational waves. In this paper, 
we consider the generic properties of gravitational waves produced during preheating, perform de- 
tailed calculations of the spectrum for several specific inflationary models, and identify problems 
that require further study. In particular, we argue that if these gravitational waves exist they 
will necessarily fall within the frequency range that is feasible for direct detection experiments - 
from laboratory through to solar system scales. We extract the gravitational wave spectrum from 
numerical simulations of preheating after X(f> 4 and micf> 2 inflation, and find that they lead to a grav- 



itational wave amplitude of around D. gw h 



10" 



This is considerably higher than the amplitude 



of the primordial gravitational waves produced during inflation. However, the typical wavelength of 
these gravitational waves is considerably shorter than LIGO scales, although in extreme cases they 
may be visible at scales accessible to the proposed BBO mission. We survey possible experimental 
approaches to detecting any gravitational wave background generated during preheating. 



I. INTRODUCTION 

In recent years, the Cosmic Microwave Background 
(CMB) has been our primary window into the primor- 
dial universe. Scale invariant density perturbations, 
sourced by quantum fluctuations during inflation and 
processed by the photon-baryon plasma, are by far the 
most satisfying explanation for the observed temperature 
anisotropies in the CMB. Inflation also predicts the exis- 
tence of a scale invariant spectrum of primordial gravita- 
tional waves, sourced by the same quantum fluctuations 
that underlie the scalar perturbations. Gravitational 
waves are only weakly coupled to matter fields, and move 
freely through the universe from the moment they are 
produced. They are thus the deepest probe of the early 
universe of which we currently know 0, H, 0, 0, H, H, Q • 

In the short term, the best hope for observing primor- 
dial gravitational waves is via their contribution to the 
B-mode of the CMB polarization. However, foreground 
weak lensing puts a fundamental limit on a B-mode signal 
sourced by gravitational waves [3, M, LLlJ ■ Consequently, 
attention has recently focussed upon direct detection ex- 
periments 0, El 13 El which, while enormously chal- 
lenging, may ultimately be more sensitive to a stochastic 
gravitational wave background than the CMB B-mode. 
These experiments are sensitive to frequencies far higher 
than those that contribute to the CMB, since the phys- 
ical sizes of detectors necessarily range from laboratory 
scales through to solar system scales. However, as the in- 
flationary spectrum is almost scale- free, the amplitude at 
short scales is not dramatically different from that seen 
at CMB scales, at least for simple models of inflation. 

The primordial spectrum is not alone: gravitational 
waves are produced whenever there are large, time- 
dependent inhomogeneities in the matter distribution. 
These more pedestrian gravitational waves are generated 
"classically" , much like radiation produced by the oscilla- 



tion of an electron, whereas the inflationary tensor per- 
turbations are sourced by quantum fluctuations in the 
spacetime background. During the evolution of the uni- 
verse, there are several hypothetical mechanisms which 
would generate large, local inhomogeneities. These in- 
clude first order phase transitions [14| , brane- world mod- 
els 0, networks of cosmic strings 0, or inho- 
mogenous neutrino diffusion |18|. Here, we investigate 
a further possibilit y: g ravitational waves produced dur- 
reheating H M Ejl ElEl El El & El El H 



Liil IMLM IM IM |35l |3l |37l IM 1391 140(, a period of 
non-thermal evolution following the end of inflation. 

While the detailed dynamics of resonance and pre- 
heating is a complicated, nonlinear problem, the ba- 
sic picture is very simple. For a given combination 
of fields, parametric resonance occurs when some sub- 
set of Fourier modes have exponentially growing solu- 
tions, driven by the oscillating inflaton field. The res- 
onant modes are quickly pumped up to a large ampli- 
tude <I> 2 w (x 2 )- Resonance ends when nonlinearities 
render further growth kinematically expensive, leaving 
the universe far from thermal equilibrium. Thcrmaliza- 
tion occurs as the excited modes dissipate their energy 
over the entire spectrum via self-interaction. The rapid 
rise in the mode amplitude can be associated with an 
exponentially growing occupation number (at least for 
bosonic species). If one transforms into position-space, 
the highly pumped modes correspond to large, time de- 
pendent inhomogeneities, ensuring the matter distribu- 
tion has a non-trivial quadrupole moment, sourcing the 
production of gravitational radiation. 

This topic was first addressed by Khlebnikov and 
Tkachev [27]], and has not been widely discussed within 
the experimental community. We believe that the time is 
ripe for revisiting this question. Since |2?| was written, 
technological approaches to gravitational wave detection 
have advanced considerably. Moreover, there have been 
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significant advances in the theoretical understanding of 
preheating in more complicated inflationary models. Fi- 
nally, numerical simulations benefit from the gains in 
computational power over this interval. It will turn out 
that for "typical" cosmological parameters, gravitational 
radiation sourced by preheating has a peak frequency in 
the MHz band. Coincidentally, a new generation of de- 
tectors has been proposed which is tuned to gravitational 
waves in this range [4j]], although their strain sensitivity 
would need improve over current values by around 10 6 
in order to detect the spectra we compute here. This is 
clearly ambitious, but probably not unreasonable when 
compared to the 25 year interval anticipated before BBO 
[Big Bang Observer] will be in a position to detect the 
inflationary gravitational wave spectrum [42j . 

Since preheating occurs in many (although not all) in- 
flationary models, any gravitational wave signal associ- 
ated with preheating would provide a new and currently 
unexplored window into inflationary physics. It is worth 
emphasizing that the physical principles that underlie the 
calculations in this paper are well established and do not 
rely on exotic physics - in particular the mechanism that 
governs the generation of the gravitational waves is sim- 
ply the usual quadrupole related emission. 1 Finally, this 
signal carries information about the epoch at the end of 
inflation, opening a new window into the early universe. 
With the stakes this high any possible observational op- 
portunity must be carefully explored. 

Unlike the primordial spectrum, gravitational waves 
induced by preheating (or any subsequent phase transi- 
tion) will not be scale-free; indeed their spectrum is in- 
dicative of the complicated processes that generate them. 
This is both a boon and a possible pitfall. A scale de- 
pendent spectrum necessarily contains more information 
than a scale-free spectrum, so the detection and subse- 
quent mapping of a gravitational wave background in- 
duced by preheating would yield a rich trove of informa- 
tion. However, gravitational wave detectors are subject 
to physical limitations, since they all ultimately measure 
deformations in (an array of) physical objects induced 
by passing gravitational waves. Thus, even in principle, 
it is hard to imagine the direct detection of gravitational 
waves below atomic scales, or beyond solar system scales. 
To be sure, this is a large range but it is very much shorter 
than any relevant cosmological scale. Fortunately, by 
a happy numerical coincidence, any gravitational wave 
spectrum generated during preheating will peak at scale 
between a few meters and millions of kilometers - which 
overlaps with the "golden window" open to direct detec- 
tion experiments. 

In Section ljTT|l . we briefly review preheating and dis- 
cuss previous work. Following that, in Section we 



1 The possibility that preheating dynamics are directly affected 
by "ba ck-reaction" from metr ic perturbations has also been dis- 
cussed |3(1 1 i ll \.'>A I i l |3S| but we do not address this effect 
in this work. 



demonstrate that preheating gravitational waves from in- 
flation scales ranging from TeV to the GUT scale will 
peak around 1 Hz to 10 8 Hz. In Section (|IV|1 we discuss 
our methodology for numerically computing the gravi- 
tational wave spectrum. In Section (|V|) we present our 
results for A</> 4 and ro 2 <ft 2 inflation, confirming and ex- 
tending the results of [23]. We note that while the in- 
flationary dynamics of the two models are very similar, 
their resonant behaviors diverge considerably. Finally, 
we discuss the theoretical and observational implications 
of our results and lay out future plans to further improve 
the computational methodology in Section IjVIfl . 



II. PREHEATING 

A key problem facing any inflationary model is to en- 
sure that inflation ends. This issue is highlighted in 
Guth's foundational paper, describing what is now known 
as old inflation, which is driven by a meta-stable false 
vacuum does not successfully terminate |43|. In new or 
chaotic (slow-roll) inflation [44L l45l I4r3 | , it was thought 
that inflation was followed by a period of reheating, where 
energy slowly bleeds from the inflaton field as it oscil- 
lates about the minimum of its potential 01 • Generating 
the correct perturbation spectrum typically requires that 
inflaton self-coupling is extremely weak, and this small 
coupling must be protected from loop corrections. Con- 
sequently, the coupling between the inflaton and other 
particles is necessarily tiny (< 10 -6 for A^ 4 ), ensuring 
that reheating proceeds slowly. 

Preheating provides a vastly more efficient mechanism 
for extracting energy from the inflaton field, and it pro- 
ceeds non-perturbatively and non-thermally via a process 
known as parametric resonance [U |3| • This is akin to 
stimulated emission in a laser: during preheating, indi- 
vidual momentum modes of fields coupled to the inflaton 
(or the inflaton itself, in some cases) have exponentially 
growing amplitudes. Consider the following action 
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As usual the Hubble parameter, H and scale factor, a, 
are related by H = a/a and during inflation the dynamics 
are described by 
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where V = V(<p) + l/2g 2 (j> 2 x 2 is a shorthand for the po- 
tential terms in For simplicity we consider only cou- 
plings to scalar fields although fermionic preheating has 
also been investigated [2^, |2g . 

In this action, <f> is the inflaton and the inflationary 
dynamics are fixed by the potential V ((/)), which is as- 
sumed to possess a minimum. At the end of inflation, 
the potential energy of the field is quickly converted into 
kinetic energy, and the field oscillates with frequency 
uifp = y/ d 2 V (4>) I d(f> 2 evaluated at the minimum of V(^). 
The solution for cj> is then approximately 



(j)(t) — sinm^t. 



(6) 



Meanwhile, the equation of motion for the \ field after- 
expanding it in Fourier modes is |24j 



Xk + 3H X 'k + (A{k) - 2qcos2z) X k = 0. 



(7) 



In the limit where we can ignore the expansion of the 
universe this is the Mathieu equation which is simply a 
harmonic oscillator with a periodic forcing function. Here 
we have made the identification A(k) — k 2 / (ni^a 2 ) + 2q, 
rescaled the time to z = m^t and used a prime to denote 
differentiation with respect to z. The crucial resonance 
parameter is 



q = ?&/(4m%) 



(8) 



The Mathieu equation possesses both oscillatory and ex- 
ponential solutions; for each individual mode k one can 
compute A(k) and q to determine whether or not it goes 
into resonance 0. Roughly speaking, for broad res- 
onance where a large number of modes are excited we 
need q > 1 HU- 

As first described by |U > some Xk will have expo- 
nentially growing solutions for realistic parameter values. 
A full treatment of parameteric resonance in an expand- 
ing universe is complicated, but we can make several 
generic statements. Firstly, preheating is very efficient 
and proceeds much more rapidly than reheating, which 
relies on tree level couplings between the inflaton and 
other matter fields (which, at minimum are provided by 
gravitational interactions). Parametric resonance typi- 
cally lasts less than a Hubble time, and in some models 
will complete in a few oscillations of the inflaton. This 
is because the resonant modes are rapidly pumped up to 
an amplitude (x) 2 ~ & 2 , cutting off resonance as it be- 
comes kinematically expensive. Once preheating ends, 
the pumped-up modes dissipate their energy via self- 
interaction with other modes, thermalizing the universe. 

Preheating leads to an initially non-thermal distribu- 
tion of energy in the Xk states. At high frequencies, cor- 
responding to Fourier modes that are much shorter than 
the size of the post- inflationary Hubble horizon, the effec- 
tive mass of the Xk state is much larger than the function 
amplitude and resonance does not occur. Meanwhile, at 
low frequencies, causality ensures that modes longer than 



the Hubble horizon are unlikely to be in resonance. 2 Con- 
sequently, we expect the spectrum to be narrow and cen- 
tered around a wavelength which is dependent on the en- 
ergy scale at the end of inflation. For the models we study 
in detail here, the gravitational wave spectrum induced 
by preheating peaks at scales 1^2 orders of magnitude 
shorter than the energy scale at the end of inflation. This 
comoving scale can be converted into a physical scale in 
the present universe once the post-inflationary behavior 
of the scale factor a(t) is specified. In Section (|lll|l we 
show that the peak wavelength has a physical wavelength 
that scales as H P 



lo oc 



1 



1 



(9) 



where the e subscript denotes the value of 4> and H at the 
end of inflation. Lowering the inflation scale reduces the 
reheating scale, and reddens the gravitational wave spec- 
trum. If the longest possible modes are excited in GUT 
scale inflation, the signal peaks around 10 7 ~ I0 8 Hz, al- 
though the excited modes are generally slightly shorter. 
As we will see below, this is a very challenging frequency 
range for any direct detection experiment. Reducing the 
inflation scale pushes the signature towards more eas- 
ily observable frequencies. Given that parametric reso- 
nance naturally cuts off at both small and large scales, 
the spectrum of any gravitational waves will cover a fixed 
range of wavelengths. As the power is thus restricted to a 
relatively narrow band, the total gravitational radiation 
remains safely below the bound from big bang nucleosyn- 
thesis m. 

The bottom line is that for a short moment the uni- 
verse is highly inhomogenous, providing a fertile ground 
for the generation of gravitational waves. Needless to say, 
preheating is a highly non-linear process and analytical 
estimates can only take us so far. Fortunately, given an 
action such as JTJ , this is a problem that can be solved nu- 
merically; we simply derive the equations of motion and 
evolve them numerically on an expanding lattice. This 
mimics the growth of the universe, but ignores the back- 
reaction of metric perturbations on the field evolution - 
an assumption that is self-consistent, as while Sp/p can 
be large, the metric perturbations are typically small. We 
use a modified version of the publicly available package 
LatticeEasy [Tjj for the numerical computations. 

To our knowledge, the generation of gravitational 
waves by preheating has been thoroughly examined only 
once before, by Khlebnikov and Tkachev j2^]. Their work 
is the starting point for this paper: we elaborate and ex- 
pand upon their treatment, considering a broader range 
of models, and taking recent developments in detector 
technology into account. In particular, preheating can 



2 It is actually not impossible to have resonance at very small k, 
but it does n ot occur in generic models of preheating. For a 
summary see l37l . 
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occur at a very broad range of scales, for example via 
hybrid inflation |2(J, |28|, |5(J] . If the scale gets low enough, 
the peak wavelength can be close to the scales probed by 
next generation observatories such as BBO 1421. In Sec- 
tion (0, we reproduce numerical results of |27J for the 
A</> 4 model, and present new results for the m 2 (f> 2 model. 
While the inflationary behavior of these two models is 
similar, their resonance structure is very different, pro- 
viding a useful crosscheck on the generality of our an- 
alytic estimates. In future work we will extend these 
numerical calculations to hybrid inflation and fermionic 
preheating. 



III. PEAK WAVELENGTH AND AMPLITUDE 

We begin our detailed analysis with a general discus- 
sion of the different parameters that determine the am- 
plitude and wavelength of any gravitational waves gen- 
erated during preheating. Consider the "usual" gravita- 
tional wave power spectrum generated by quantum fluc- 
tuations of the background, 



mode k ~ H r 



where V^Jj is the energy scale of inflation and fl r h 2 w 
4 x 10~ 5 is the total density of radiation today. The 
effective number of degrees of freedom in the radiation 
at matter-radiation equality and today are given by 
and go respectively. This form of the power spectrum 
is slightly non-standard as tensor modes are usually ex- 
pressed via Ph = 8ttH 2 /m| r There are two salient fea- 
tures to this spectrum. The first is that it is (almost) 
scale-invariant, since each mode is frozen out at an ap- 
proximately constant energy scale, namely the inflation 
scale. The second is that the power is minute; for the 
most optimistic scenario where inflation occurs around 
the GUT scale, fl gw ,i n f(k)h 2 < 10~ 14 . The current up- 
per limit on the scale of inflation from WMAP observa- 
tions for single-field inflation models is U 1//4 < 3.3 x 10 16 
GeV, corresponding to fl gw s n fh 2 < 2 x 10" 15 [TM l5lj| 

The gravitational waves which we are considering in 
this paper are not directly sourced by quantum fluctua- 
tions; instead they are generated by the classical motion 
of particles during preheating. As is well-known, accel- 
erated motion generates a quadrupole moment, leading 
to the generation of gravitational radiation. During pre- 
heating at the end of inflation, large inhomogeneities in 
the matter fields are generated by the selective pump- 
ing of modes in parametric resonance. These large in- 
homogeneities, as first shown in 27] for the A</> 4 model, 
are sufficiently large to produce gravitational waves with 
amplitudes many orders of magnitude larger than those 
produced by the quantum fluctuations. An analytical es- 
timate |27^ . for an inflaton with an effective oscillation 
frequency m, coupled to a massless scalar field \ with a 
g 2 (j> 2 X 2 term, yields the peak amplitude at the resonance 



32 / V. 



1/3 



(10) 



n gw (k ~ H e )h 2 w n r h 2 



g 2 m 2 l 



1/3 



(11) 



In other words, the amplitude probes the oscillation scale, 
in contrast to the primordial spectrum which probes the 
inflation scale. If we plug in the usual field values for 
chaotic inflation fh = at the end of inflation <fi « m P1 
such that V e = m 2 </> 2 /2 = m 2 ^m 2 , i /2, we see from equa- 
tion (|10|) that the amplitude of the gravitational waves 
generated by preheating is 1 / g 2 larger than the inflation- 
ary spectrum. This is a significant boost, as we expect 
g 2 < 1CU 6 . Note that in models such as hybrid inflation 
3 one can decouple the oscillation frequency m from the 
inflaton mass m, and this simple relationship has to be 
revisited polEs)]. 

However, only a finite range of modes excited during 
preheating. If the power was generated at scales corre- 
sponding to, say, atomic distances today, then our hope 
of detecting any gravitational waves induced by preheat- 
ing would be dashed. On causal grounds, we expect that 
resonant modes have a wavelength roughly equal to or 
less than the Hubble length at the end of inflation, 1/H e : 



(12) 



where V e is the inflationary potential. After inflation, 
the universe reheats to a temperature T rh . During the 
subsequent radiation dominated phase, the Hubble pa- 
rameter scales as H — iJ»(a*/a e ) 2 until matter-radiation 
equality at T*. Meanwhile, the scale factor evolves as 
a* = ao(go/<?*) 1 ^ 2 (2o/T'*) from matter-radiation equal- 
it y u ntil today when a = 1. Thus for a physical length 
I |27| and physical wavevector k we have 



1 1/2 

1 9* 

fc 1/3 

So 



8tt 3 
~9CT 



-1/4 



y/H e Mp 



0.5 



f = 6 X 10 



10 



Hz 



(13) 



(14) 



where we have used go/g* — 1/100 in the second line. 
Plugging in the lowest excitable frequency k — H e , where 
we more or less expect peak gravitational waves pro- 
duction to occur, we obtain the scaling relation J^J, as 
claimed earlier. 



3 We note that in hybrid inflation, preheating amplification of the 
perturbations is achieved through a combination of parametric 
resonance and tachyionic instabilities. We thank Gary Felder for 
pointing this out to us. 



5 



The inverse scaling is particularly important: it means 
that the pertinent wavelengths are longer for smaller in- 
flationary energy scales. If we assume instantaneous re- 
heating after inflation for GUT scale inflation H e as 10 13 
GeV, I w 1 - 10 meters, and / w 10 7 to 10 s Hz. Lower- 
ing the inflationary scale reduces power in the primordial 
gravitational wave spectrum making it harder to detect, 
as quantified by equation (jTU|l . However, this also red- 
dens the peak power of any preheating generated gravita- 
tional waves, making them easier for us to observe. This 
follows because the strain sensitivity ht of a detector 
scales as fl gw /f 3 j4(|, i.e. for the same value of £l gw we 
have to build a more sensitive detector if the frequencies 
are higher. In addition, if inflation occurs at a lower scale, 
then the gravitational wave energy density will be diluted 
less by expansion following preheating, again increasing 
our chance of observing them. On the other hand, if the 
gravitational waves are generated at a lower scale the off- 
diagonal terms olT^ v are smaller for fixed 5p/p, and will 
be less efficient sources of gravitational radiation. On the 
basis of the limited calculations performed in this paper, 
we see that the last two effects roughly cancel and D, gw 
does not depend strongly on H e . However, further work 
will be needed before we can safely say that this is true 
of all models which undergo preheating. 

In this naive analysis, a few subtle points have been 
glossed over. The peak resonance modes are usually not 
exactly at the Hubble scale; instead they are frequently 
1-2 orders of magnitude smaller [2l[ |H |2<|. This 
has the effect of pushing the observable modes to a bluer 
band. On the other hand, preheating does not always 
start immediately after inflation ends; peak particle pro- 
duction occurs when the amplitude of the field perturba- 
tions 5<fi/4> grows to order unity, which need not happen 
quickly. In the models looked at here, the Hubble param- 
eter during peak gravitational wave production is about 
1 — 2 orders of magnitude smaller than H e , shifting the 
observable modes to a redder band. 

It is also worth noting that the gravitational waves are 
generated causally within the Hubble volume, and thus 
the phases of the individual modes are uncorrelated - 
unlike the primordial spectrum. This is a generic feature 
of all causally generat ed p erturbation spectra, and is a 
powerful discriminant |52j. Unfortunately, direct detec- 
tion experiments cannot dinstinguish the coherence (or 
lack of) of the gravitational waves as their signal is an in- 
tegral over some time interval greater than the frequency 
scale. To do this, one must find a processed imprint on 
a fixed time-slice. 

While there is an upper bound on the inflationary en- 
ergy scale from the contribution of tensor modes to the 
CMB, the lower bound is very weak. At minimum, the 
post-inflationary universe must be hot enough to per- 
mit baryogenesis and nucleosynthesis. We conservatively 
assume that the former occurs via electroweak scale pro- 
cesses, so we can easily have Ve^ 4 as low as the TeV 
scale. Nuclear reactions necessarily take place at MeV 
scales and ensuring successful nucleosynthesis provides 



an absolute lower limit on the reheating temperature. 
This corresponds to gravitational wave peak wavelength 
scales ranging from laboratory scales through to solar 
system scales. 

IV. GRAVITATIONAL WAVE PRODUCTION 

Equation (jTT|l suggests that the preheating induced 
gravitational wave spectrum is larger than its primor- 
dial counterpart. To obtain an actual power spectrum, 
the highly nonlinear physics of preheating forces us to 
turn to numerical methods. We use LatticeEasy 
to simulate the evolution of the early universe, solving 
the equations of motion for a set of interacting scalar 
fields in a flat Friedmann-Robertson- Walker (FRW) Uni- 
verse. The fields become highly inhomogeneous, which 
is important for the generation of gravitational waves. 
This does not immediately make the metric perturba- 
tions large [53|. Consequently, we can solve the non- 
linear field evolution numerically while assuming a rigid 
spacetime background, and then extract the spectrum of 
gravitational radiation produced during preheating. 

We extended LatticeEasy to compute the gravita- 
tional wave spectrum generated duringpreheating. We 
follow the approach of j23| (see also [3]), reproducing 
their results for the quartic A</> 4 /4 + <7 2 2 y 2 /2 model. In 
addition, we compute the gravitational wave spectrum 
for the quadratic inflation model m 2 </> 2 /2 + g 2 </> 2 x 2 /2. 
We leave the simulation of other models such as the neg- 
ative coupling — g 2 2 x 2 /2 |2l| or hybrid inflation models 
to future work. 

We now sketch the approach we use to compute the 
spectrum. We begin by considering the energy radiated 
in gravity waves in a frequency interval diu and a solid 
angle dfl, given by 

^ = 2GK l] . lm uj 2 T^*{^u J )T lm {^u)dw (15) 

where T 4J ' is the stress tensor describing the source mat- 
ter fields. Here i,j run over the spatial indices and the 
projection tensor is given by |54j 

A-ij,lm{ty — fiij^lm 2kjk rn 5u + —kikjklkm 

— \hj5lm + -^Sijklkm + -Sjihkm (16) 

with unit vector k = k/w. Strictly speaking, this for- 
mula is only valid for linearized gravity in Minkowski 
space. A more accurate calculation will involve solving 
the equations of motion for linearized gravity on a curved 
background. 

To see why the use of this formula is justified, consider 
the gravitational wave energy emitted by a three dimen- 
sional box of conformally flat spacetime with physical 
size I x I x I. The energy density in this box is 

dp{uj) = 87rGr 3 Ay, im w 2 T^'*(|k|,u;)T' m (|k|,w)dw (17) 
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where we have assumed that the spectrum is isotropic. 4 
From causal arguments alone, only modes of wavelengths 
equal to or shorter than 1 / H will be generated, imposing 
a natural cut-off at long scales. Thus, provided we choose 
I > 1/H we will effectively capture the essential physics. 
Depending on how efficient preheating is in a particular 
model, the entire phase can last for several Hubble times. 
However, the gravitational waves are produced near the 
end of preheating, as the inhomogeneities in the fields be- 
come large. We therefore expect the gravitational wave 
power to be generated in a short burst, and numerical 
simulations confirm this suspicion. Thus it is a good ap- 
proximation to assume that the gravitational wave source 
is localized in the box |55l ] . 

In our simulations, we begin our computations at the 
end of inflation, near the beginning of the parameteric 
resonance phase. We end our simulations when the fields 
are stabilized and parameteric resonance ends. We subdi- 
vide the the spacetime into discrete 4-D boxes of spatial 
sizes L 3 and time interval r = L, where L and r are 
the conformal length and time respectively. Our phys- 
ical box size thus scales roughly as a 3 . The choice of 
t = L is purely operational, allowing us to fix our Fourier 
variables to be the conformal frequency and conformal 
wavevector such that |k co „/| = u> con f. One can in prin- 
cipal decompose the conformal time differently, but that 
would unnecessarily complicate matters. We fix L so 
that during the period of gravitational wave production 
ah > 1 /H and the box is larger than the effective Hubble 
horizon. 

We assume that each "box", labeled a, is a localized 
source, and compute the total gravitational wave density 
produced for each box using equation l|17J) . We then 
sum them up, diluting them appropriately as follows 

dp gw {a e ) _y^ dpgw(a a ) ( Q e \ 4 
d In u) d\nuj \a a J 

a x / 

where a a is the scale factor taken at the middle of the 
box in conformal time and a e is the scale factor at the 
end of inflation. Meanwhile, for each box a 

%M = ^Gu% 3 K, hlm T^T lm (19) 
din a; 

where l 3 ^ is the physical size of the box at time a a and 
lu is the physical frequency. 

Finally, putting everything together, the total density 
of gravitational waves today is given by 

^ 2 =^ 2 %^f^y /3 (20) 

a In ui \g*J 



4 Isotropy allows us to choose any direction for the projection vec- 
tor Ay i m : we chose for simplicity k = (1, 0, 0). We checked that 
this assumption is robust by showing that the final simulation 
results are not sensitive to different choices of direction. 



Log[Q g „ h 2 ] 




FIG. 1: The gravitational spectrum for the A<^ 4 model with 
A = 1(T 14 and g 2 /\ = 1.2 (dash line) and 120 (full line) 
respectively . As expected, it is peaked around 10 7 ~ 10 8 Hz 
and spans about 2 decades. The horizon size at the time of 
preheating imposes the low frequency cut-off, while the high 
frequency cut-off is due to the fact that high momentum \ 
particles are energetically too expensive to be created. Notice 
that the power is roughly inversely proportional to g 2 . 

We should mention that in equation IjlSJ) and hence 
equation JSUJ), we have implicitly assumed that the uni- 
verse is radiation dominated at the end of preheating, 
which is not true for certain chaotic models. 



V. NUMERICAL RESULTS 

In this section, we give numerical results for the gravi- 
tational wave spectrum produced during resonance in two 
different models: A0 4 and m 2 4> 2 - While the inflationary 
dynamics of these two systems are very similar, there 
is considerable divergence in the resonance structure be- 
tween the models, making this a useful generalization of 

m 

A. Quartic Inflation (A0 4 ) 

To test our code, we reproduce Khlebnikov and 
Tkachev's results 27] for \<j) 4 with a 4> 2 x 2 term 

V(0,X) = ^ 4 + ^W- (21) 

From the perspective of preheating, this model is atypical 
[22j as it possess only a weak resonance band. Even so, 
we still sec significant production of gravity waves. 

Following [231, we set A = 10~ 14 and g 2 /X = 120, cor- 
responding to a resonance parameter q w 120 from equa- 
tion (JHJ. In this model, inflation ends around the GUT 
scale, where ^>o ~ m Pi, ° 1 ' H en d w 10 12 GeV. We begin 
our simulation on a 256 3 size lattice from that time and 
run it until preheating ends around H ss 10 7 GeV. Pa- 
rameteric resonance peaks around H vea k ~ 10 8 GeV, and 
the size of the box is chosen to ensure that its physical 
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variance 
r 



variance 
r 





FIG. 2: The variance of the cj> (dash) and \ (Ml) fields for a 
A</> 4 inflation model with g 2 /\ = 120. The time coordinate is 
the conformal time r in units of 2.5 x 10 _12 GeV _1 while the 
variances are in m 2 , units. The associated hubble parameter 
during the rapid rise while in major resonance phase is H rts 
10 GeV, thus this phase lasts less than a Hubble time. 



size at this time I as 1/H pea k- With a A</> 4 potential, the 
background spacetime scales like a radiation dominated 
universe during parametric resonance. 

Using i|14|) . the present frequency associated with the 
Hubble parameter during preheating is 10 6 Hz. From 
figure OJj we see that the peak frequency is actually 
10 7 ~ 10 8 Hz, suggesting that the peak resonance modes 
are about two orders of magnitude smaller than the Hub- 
ble wavelength. The amplitude of the gravitational waves 
peaks at around fl gw « 10~ 9 , consistent with equation 

03. 

From the plot, we see that even at the lower end of the 
relevant frequency range which is easier to detect, Q gw ~ 
10 -11 . This is 3 orders of magnitude larger than the 
primordial spectrum. Beyond that at lower frequencies, 
we expect the spectrum to undergo a steep A; 3 decline. 
This k 3 superhorizon tail is a common property for a 
spectrum which is causally generated inside the Hubble 
horizon |34j . In cases where the inflationary scale and 
thus the intrinsic stochastic background is very low (so 
it does not mask the signal) this fc 3 tail might be easier 
to detect than the peak wavelengths, given the physical 
limitations on realistic detectors. 



B. Quadratic Inflation (mi<^ 2 ) 



FIG. 3: The variance of the 



inflation model with q - 



We now turn to the m 2 ^ 2 model |4fj 
interaction term as before 



with the same 



V(0,x) = im^ 2 + I ff V x 2 - 



(22) 



The amplitude of the CMB temperature anisotropy re- 
quires rricj, as 10 13 GeV. At the end of inflation <j> as m Pl . 
Choosing g 2 — 2.5 x 10~ 7 , gives a resonance parameter of 
q as 2.5 x 10 5 , via JSJ. Again we begin our simulation af- 
ter inflation ends at H en d as 10 13 GeV through the peak 



is the cosmic time t in units of m^ 1 



■ (dash) and x (full) fields for 
2.5 x I0 5 . The time coordinate 
S.33 x lO-^GeV- 1 

V 

while the variances are in m% l units. The associated hub- 
ble parameter during the rapid rise while in major resonance 
phase is H as 10 12 GeV, thus this phase lasts much less than 
a Hubble time. 



Log[Q g „ h'J 
-10 r 

-10.5 

-11 ■ 
-11.5 

-12 ■ 
-12 . 5 

-13 ■ 




Log (f/Hz) 



FIG. 4: The gravitational spectrum for a m|0 2 inflation 
model, with parameters g 2 — 2.5 x 10~ 7 and = I0~ 6 m P i ~ 
10 13 GeV. The resonance parameter here q — 2.5 x I0 5 . The 
slight rise at high frequencies after the sharp drop is a numer- 
ical artifact. 



preheating phase at H pea k « 10 11 GeV, until the end of 
preheating. Since is non-zero, the universe evolves as 
if it was matter dominated 5 with the scale factor grow- 
ing 30-fold. At the end of preheating, we assume that the 
universe reheats normally and enters a radiation domi- 
nated phase. This is in principle not a valid assumption, 
as some numerical results have shown that it is difficult 
for quadratic inflation to reheat to radiation domination 
without a trilinear coupling [HiJ. In this paper, we are 
using it as a toy model to illustrate preheating for a in- 
flation model with a different mass scale. 

We simulated this model on a 256 3 lattice, with the 



The equation of state w = (p) / (p) fluctuates rapidly between 1 
and —1 around a center value of 0. 
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Log[Q g .„ h 2 ] 
-10 r 




FIG. 5: The gravitational spectrum for a rn^tj) 2 inflation 
model, with parameters g 2 = 2.5 x 10 -7 and ~ 10~ 7 m P i = 
10 12 GeV. The resonance parameter here q = 2.5 x 10 5 . Al- 
though this lower mass model is ruled out by the CMB, it 
serves as a useful prototype to illustrate that the peak will be 
reddened for a lower scale inflation. 

results as shown in figure (J2J. 

Using (|llf> . we expect ft gw h 2 ss 10~ 10 at peak which 
matches the result of our detailed calculation. Although 
inflation ends at H en d ~ 10 13 GeV, peak resonance oc- 
curs at H pea k ~ 10 11 GeV, which sets the comoving size 
of our lattice. Figure (@J shows that the peak frequency 
is actually 10 8 ~ 10 9 Hz, a couple of orders of magnitude 
smaller than the Hubble parameter during reheating. 

Finally, we present the results for a model with a lower 
mass, Tn^ = 10 12 GeV in figure 0. This model is ruled 
out by CMB data, but it demonstrates the way in which 
the gravitational wave spectrum generated by preheating 
depends on the inflationary scale. In this model H « 10 12 
GeV, and as expected from equation the peak loca- 
tion is reddened by a factor of y/\Q. The observable power 
remains comparable to the previous model. The emitted 
power is reduced, but since the overall expansion of the 
universe is reduced by the lower reheating temperature, 
the values of Q gw h 2 today is roughly fixed. It is tempt- 
ing to conjecture that the cancellation between these two 
effects will be seen in other preheating induced gravita- 
tional wave spectra, and that the VLg W h 2 ~ 10~ 10 seen 
here will prove to be a generic value j^ll E3| . 



VI. SUMMARY AND FUTURE PROSPECTS 

We have carefully investigated the production of grav- 
itational waves during prehea ting , reproducing the work 
of Khlebnikov and Tkachev (23 for the A</> 4 inflation 
model and extending it to the m 2 2 case. For both 
models we show numerically that preheating is a sizable 
source of gravitational waves with frequencies of around 
10 6 - 10 s Hz, and peak power of Q gw h 2 » lfr 9 - 1CT 11 . 
We present simple scaling arguments to predict the over- 
all properties of the spectrum for a broader class of in- 
flationary models. We see that the spectrum of grav- 



itational waves induced by preheating peaks at a scale 
proportional to 1/yH, where H is the Hubble param- 
eter during preheating, and generally somewhat smaller 
than the scale of inflation. Thus, lowering the inflation- 
ary scale reddens the spectrum and makes it easier to 
observe. This is in contrast to the primordial inflationary 
spectrum, which is roughly scale invariant and becomes 
harder to observe as inflationary scale is lowered. 

We now ask what we can learn about inflation if we 
detect a spectrum of gravitational waves generated dur- 
ing preheating. Its two most basic features, the peak 
frequency and the amplitude, represent the reheat and 
the oscillation scales respectively. As the reheat scale is 
lower than the inflation scale, its detection would impose 
a lower bound on the inflationary scale. Its usefulness 
as a probe of inflation is amplified if we have a separate 
probe of the scale of inflation, say from the CMB B-mode 
observations. 

The oscillation scale is harder to interpret, as it is of- 
ten highly model dependent. For single scalar field infla- 
tion, such as the models considered here, knowledge of 
the reheat scale would constrain the coupling constant 
g 2 . More optimistically, the structure of the spectrum 
encodes information about resonance and preheating, so 
if we can predict the structure accurately we potentially 
probe the detailed mechanics of preheating. Such an en- 
deavour will require more careful computations and sim- 
ulations than we present in this paper. 

Further progress on this problem can be made in two 
ways The first is to further refine the code to ac- 

commadate a larger class of models, particularly hybrid 
inflation models which have an essentially arbitrary infla- 
tionary scale. Alternatively, as alluded to in Section 1(1 V^. 
a more sophisticated theoretical calculation would be to 
directly solve the evolution equations for the off-diagonal 
parts of the perturbed Einstein tensor, which are sourced 
by Tij. This would avoid any ambiguity concerning the 
use of a formula that is only strictly applicable in flat 
space, and it would avoid the need to run the code for 
a finite number of "boxes" , since we would only need to 
take a Fourier transform at the end of the computation. 

Investigating these gravitational waves is timely, since 
there is currently considerable interest in the direct detec- 
tion of gravitational waves. At the moment, several pro- 
posals based on different technologies are being actively 
pursued. At the solar-system scale, the space-based inter- 
ferometer LISA (53 will probe frequencies from 10~ 2 Hz, 
which is probably too small for the gravitational waves 
we are considering. The proposed BBO \^ and also 
the Deci-hertz Interferometer Gravitational Wave obser- 
vatory (DECIGO) 0| missions are sensitive to frequen- 
cies on the order of 1 Hz, and would probe gravitational 
waves arising from preheating after TeV scale inflation. 

The array of terrestrial interferometers also probes fre- 
quencies corresponding to preheating following low-scale 
inflation. These experiments include GEO600 60], LIGO 
[HI, TAMA IH and VIRGO [63|. These are sensitive 
to scales between 100 ~ 1000 Hz and may be able to 
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FIG. 6: Projected sensitivities of gravitational wave detectors (dotted lines) to a stochastic gravitational wave spectrum (full 
lines) in log(S7 sm /i 2 ) and log(//Hz). The limits are for correlated observatories, except for LISA. The primordial inflation 
spectrum is plotted for a scale-invariant GUT scale inflation. The gravitational wave spectra generated by the quartic and 
quadratic type inflation models for a GUT scale inflation are schematically shown here; see Section for more detailed plots. 
In addition, we have extrapolated a f 3 tail to the red of the simulated results, which is a conservative estimate on the degree of 
the Hubble scale cut-off to the gravitational wave spectrum [34|. We have also added the expected peak point of a TEV scale 
hybrid inflation model with high oscillation scale, calculated using the naive formula 111! 12q| . 



probe a stochastic background in the interesting range 
Qgwh 2 ~ 10 -10 , if they are correlated. At even higher 
frequencies in the KHz ran gje, w e have a slew of resonant 
bars detectors 0, |6^, Ir3rj l67|. Once correlated, these 
resonant bars have a potential to reach fl gw h 2 as 10~ 5 
|49| , which puts the signals we see here out of their reach. 
However, theoretical studies suggest that correlating hol- 
low spherical detectors may eventually allow us to reach 

il gw h 2 « 10" 9 m. 

At even higher frequencies from 10 3 ~ 10 5 Hz, there 
has been a proposal to build a superconducting resonant 
cavity detector called the Microwave Apparatus for Grav- 
itational Wave Observation (MAGO) [Hi! E3 Al- 
though the strain sensitivity ht ~ 10 _21 Hz -1 / 2 for the 
prototype is expected to be comparable to resonant bar 
detectors at 4 x 10 3 Hz, the large f 3 suppression from 
the relation Q gw h 2 oc hjf 3 means that at these frequen- 
cies we can only reach ft gw h 2 ss O(l). An improvement 
of 5 ~ 6 orders of magnitude in the strain sensitivity 
is needed to reach Vt gw h 2 sa 1CU 10 , a possibility which 
may be achieved by further refinements to the prototype 



and/or construction of an array of such detectors |4l| . 
By way of comparison, we note that the best hope for 
observing a primordial gravitational wave background is 
currently provided by BBO, which has a lead time of 
at least 20-25 years. In this context hoping for a large 
extrapolation of detector technologies at high frequen- 
cies is perhaps not excessively optimistic. In Figure El 
we sketch the sensitivities of the leading interferometric 
detectors along with the expected stochastic background 
produced during preheating for the models we discuss in 
detail here. 

Finally, if one was to ever make a concerted attempt 
to detect a gravitational wave spectrum associated with 
preheating, one would need to be understand other po- 
tential sources that could supply a stochastic background 
of gravitational waves. This includes any first order phase 
transition in the early universe (such as the electroweak 
scale), or decays from cosmic strings. In addition, the 
presence of a rising component in the spectrum illustrates 
the dangers of using a (locally) positive spectral index of 
any detected stochastic gravitational wave background 
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to rule out inflation in favour of alternative cosmogenesis 
ideas such as ekpyrosis 0, |?2 or pre-big bang scenarios 
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